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ABSTRACT 

We present near-IR spectroscopy of red supergiant (RSG) stars in NGC 6822, obtained with the 
new VLT-KMOS instrument. From comparisons with model spectra in the J-band we determine the 
metallicity of 11 RSGs, finding a mean value of [Z] = —0.52 ±0.21 which agrees well with previous 
abundance studies of young stars and HII regions. We also find an indication for a low-significance 
abundance gradient within the central 1 kpc. We compare our results with those derived from older 
stellar populations and investigate the difference using a simple chemical evolution model. By com¬ 
paring the physical properties determined for RSGs in NGC 6822 with those derived using the same 
technique in the Galaxy and the Magellanic Clouds, we show that there appears to be no signifi¬ 
cant temperature variation of RSGs with respect to metallicity, in contrast with recent evolutionary 
models. 

Subject headings: Galaxies: individual: NGC 6822 - stars: abundances - stars: supergiants 


1. INTRODUCTION 

A promising new method to directly probe chemical 
abundances in external galaxies is with J-band spec¬ 
troscopy of red supergiant (RSG) stars. With their peak 
flux at ^1 /im and luminosities in excess of 10 1 2 3 4 5 L©, RSGs 
are extremely bright in the near-IR, making them poten¬ 
tially useful tracers of the chemical abundances of star¬ 
forming galaxies out to large distances. To realise this 
goal, Davies et al. (2010) outlined a technique to de¬ 
rive metallicities of RSGs at moderate spectral resolv¬ 
ing power (R ^3000). This technique has recently been 
refined using observations of RSGs in the Magellanic 
Clouds (Davies et al. 2015) and Perseus OB-1 (Gazak 
et al. 2014a). Using absorption lines in the J-band from 
iron, silicon and titanium, one can estimate metallicity 
([Z] = log Z/Z©) as well as other stellar parameters (effec¬ 
tive temperature, surface gravity and microturbulence) 
by fitting synthetic spectra to the observations. Owing 
to their intrinsic brightness, RSGs are ideal candidates 
for studies of extragalactic environments in the near-IR. 

To make full use of the potential of RSGs for this sci¬ 
ence, multi-object spectrographs operating in the near- 
IR on 8-m class telescopes are essential. These instru¬ 
ments allow us to observe a large sample of RSGs in a 
given galaxy, at a wavelength where RSGs are brightest. 
In this context, the iGband Multi-Object Spectrograph 
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(KMOS; Sharpies et al. 2013) at the Very Large Tele¬ 
scope (VLT), Chile, is a powerful facility. KMOS will 
enable determination of stellar abundances for RSGs out 
to distances of ^10 Mpc. Further ahead, a near-IR multi¬ 
object spectrograph on a 40-m class telescope, combined 
with the excellent image quality from adaptive optics, 
will enable abundance estimates for individual stars in 
galaxies out to tens of Mpc, a significant volume of the 
local universe containing entire galaxy clusters (Evans et 
al. 2011). 

Here we present KMOS observations of RSGs in 
the dwarf irregular galaxy NGC 6822, at a distance of 
^0.46Mpc (McConnachie 2012, and references therein). 
Chemical abundances have been determined for its old 
stellar population (e.g. Tolstoy et al. 2001; Kirby et 
al. 2013), but knowledge of its recent chemical evo¬ 
lution and present-day abundances is somewhat lim¬ 
ited. Observations of two A-type supergiants by 
Venn et al. (2001) provided a first estimate of stellar 
abundances, finding log (Fe/H) + 12 = 7.01 ±0.22 and 
log(0/H)±12 = 8.36 ±0.19, based on line-formation 
calculations for these elements assuming local thermo¬ 
dynamic equilibrium (LTE). A detailed non-LTE study 
for one of these objects confirmed the results finding 
6.96 ±0.09 for iron and 8.30 ±0.02 for oxygen (Przy- 
billa 2002). Compared with solar values of 7.50 and 
8.69, respectively (Asplund et al. 2009), this indicates 
abundances that are approximately one-third solar in 
NGC 6822. A study of oxygen abundances in HII regions 
(Lee et al. 2006) found a value of 8.11 ± 0.10, confirming 
the low metallicity. 

NGC 6822 is a relatively isolated Local Group galaxy, 
which does not seem to be associated with either M31 
or the Milky Way. It appears to have a large extended 
stellar halo (Letarte et al. 2002; Hwang et al. 2014) as 
well as an extended HI disk containing tidal arms and a 
possible HI companion (de Blok & Walter 2000). The HI 
disk is orientated perpendicular to the distribution of old 
halo stars and has an associated population of blue stars 
(de Blok & Walter 2003; Komiyama et al. 2003). This 
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led Demers et al. (2006) to label the system as a ‘polar 
ring galaxy’. A population of remote star clusters aligned 
with the elongated old stellar halo have been discovered 
(Hwang et al. 2011; Huxor et al. 2013). In summary, the 
extended structures of NGC6822 suggest some form of 
recent interaction. 

In addition, there is evidence for a relatively constant 
star-formation history within the central 5 kpc (Weisz et 
al. 2014) with multiple stellar populations (Battinelli et 
al. 2006; Sibbons et al. 2012). This includes evidence for 
recent star formation in the form of a known population 
of massive stars, as well as a number of HII regions (Venn 
et al. 2001; de Blok & Walter 2006; Hernandez-Martfnez 
et al. 2009; Levesque & Massey 2012). 

In this paper we present near-IR KMOS spectroscopy 
of RSGs in NGC 6822 to investigate their chemical abun¬ 
dances. In Section 2 we describe the observations; 
Section 3 describes the data reduction and Section 4 
details the derived stellar parameters and investigates 
the spatial distribution of the estimated metallicities in 
NGC 6822. In Section 5 we discuss our results and Sec¬ 
tion 6 concludes the paper. 

2. OBSERVATIONS 
2.1. Target Selection 

Our targets were selected from optical photometry 
(Massey et al. 2007), combined with near-IR ( JHK S ) 
photometry (for details see Sibbons et al. 2012) from the 
Wide-Field Camera (WFCAM) on the United Kingdom 
Infra-Red Telescope (UKIRT). The two catalogues were 
cross-matched and only sources classified as stellar in the 
photometry for all filters were considered. 

Our spectroscopic targets were selected principally 
based on their optical colours, as defined by Massey 
(1998) and Levesque & Massey (2012). Figure 1 
shows cross-matched stars, with the dividing line at 
(B — V) = 1.25 x (V — R) + 0.45. All stars redder than 
this line, above a given magnitude threshold and with 
V — R > 0.6 are potential RSGs. 

Distinguishing between RSGs and the most luminous 
stars on the asymptotic giant branch (AGB) is difficult 
owing to their similar temperatures and overlapping lu¬ 
minosities. Near-IR photometry can help to delineate 
these populations, with RSGs located in a relatively well- 
defined region in the near-IR (J — K) colour-magnitude 
diagram (CMD), as discussed by Nikolaev & Weinberg 
(2000). The near-IR CMD from the WFCAM data is 
shown in Figure 2 and was used to further inform our 
target selection. Employing the updated CMD criteria 
from Cioni et al. (2014) - modified for the distance and 
reddening to NGC 6822 - all of our potential targets from 
the combined optical and near-IR criteria are (notion- 
ally) RSGs. Near-IR photometry for our targets was also 
given by Whitelock et al. (2013) from observations with 
the Infrared Survey Facility; we note that, from their 
analysis of the multi-epoch data, none of our targets was 
identified as an AGB star. 

The combined selection methods yielded 58 candidate 
RSGs, from which 18 stars were observed with KMOS, as 
shown in Figure 3. The selection of the final targets was 
defined by the KMOS arm allocation software karma 
(Wegner & Muschielok 2008), where the field centre was 
selected to maximise the number of allocated arms, with 



Figure 1. Two-colour diagram for stars with good detections 
in the optical and near-IR photometry in NGC 6822. The black 
dashed line marks the selection criteria using optical colours, as 
defined by Levesque & Massey (2012). Red circles mark all stars 
which satisfied our selection criteria. Large blue stars denote tar¬ 
gets observed with KMOS. The solid black line marks the fore¬ 
ground reddening vector for E(B — V) = 0.22 (Schlegel et al. 1998). 
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Figure 2. Near-IR CMD for stars classified as stellar sources in 
the optical and near-IR catalogues, plotted using the same symbols 
as Figure 1. This CMD is used to supplement the optical selection. 
The solid black line marks the foreground reddening vector for 
E(B — V) = 0.22 (Schlegel et al. 1998). 

priority given to the brightest targets. Optical spec¬ 
troscopy of eight of our observed stars, confirming them 
as RSGs, was presented by Levesque & Massey (2012). 

2.2. KMOS Observations 

The observations were obtained as part of the KMOS 
Science Verification programme on 30 June 2013 (PI: 
Evans, 60.A-9452(A)), with a total exposure time 
of 2400 s (comprising 8 x 300s detector integrations). 
KMOS has 24 deployable integral-field units (IFUs) each 
of which covers an area of 2 / /8 x 2 '!8 within a 7.'2 field of 
view. The 24 IFUs are split into three groups of eight, 
with the light from each group relayed to different spec¬ 
trographs. 

Offset sky frames (Of5 to the east) were interleaved 
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Figure 3. Spatial extent of the KMOS targets over a Digital 
Sky Survey (DSS) image of NGC6822. Blue filled circles indicate 
the locations of the observed RSGs. Red open circles indicate the 
positions of RSG candidates selected using our photometric criteria 
(see Section 2.1). 

between the science observations in an object (O), 
sky (S) sequence of: O, S, O, O. The observations were 
performed with the YJ grating (giving coverage from 
1.02-1.36/im); estimates of the mean delivered resolv¬ 
ing power for each spectrograph (obtained from the 
KMOS/esorex pipeline for two arc lines) are listed in 
Table 1. 

In addition to the science observations, a standard set 
of KMOS calibration frames was obtained consisting of 
dark, flat and arc-lamp calibrations (with flats and arcs 
taken at six different rotator angles). A telluric standard 
star was observed with the arms configured in the science 
positions, i.e. using the KMOSspec_caLstdstarscipatt 
template in which the standard star is observed sequen¬ 
tially through all IFUs. The observed standard was 
HIP97618, with a spectral type of B6 III (Houk & Smith- 
Moore 1988). 

A summary of the observed targets is given in Table 2. 
A signal-to-noise (S/N) ratio of > 100 per resolution el¬ 
ement is required for satisfactory results from this anal¬ 
ysis method (see Gazak et al. 2014a). We estimated the 
S/N ratio of the spectra by comparing the counts in the 
brightest spatial pixels (within the 1.15-1.22/im region) 
of each source with the counts in equivalent spatial pix¬ 
els in the corresponding sky exposures (between the sky 
lines). The S/N ratio estimated is knowingly an under¬ 
estimate of the true S/N ratio achieved. 

3. DATA REDUCTION 

The observations were reduced using the recipes pro¬ 
vided by the Software Package for Astronomical Reduc¬ 
tion with KMOS (SPARK; Davies, R. et al. 2013). The 
standard KMOS/esorex routines were used to calibrate 
and reconstruct the science and standard-star data cubes 
as outlined by Davies, R. et al. (2013). Sky subtraction 
was performed using the standard KMOS recipes and tel¬ 
luric correction was performed using two different strate¬ 


gies. Throughout the following analysis all spectra have 
been extracted from their respective data cubes using a 
consistent method (i.e. the optimal extractions within 
the pipeline). 

3.1. Three-arm vs. 24-arm Telluric Correction 

The default template for telluric observations with 
KMOS is to observe a standard star in one IFU in each 
of the three spectrographs. However, there is an alterna¬ 
tive template which allows users to observe a standard 
star in each of the 24 IFUs. This strategy should provide 
an optimum telluric correction for the KMOS IFUs but 
reduces observing efficiency. 

A comparison between the two methods in the U-band 
was given by Davies, R. et al. (2013), who concluded that 
using the more efficient three-arm method was suitable 
for most science purposes. However, an equivalent analy¬ 
sis in the Y J-band was not available. To determine if the 
more rigorous telluric approach is required for our anal¬ 
ysis, we observed a telluric standard star (HIP97618) in 
each of the 24 IFUs. This gave us the data to investigate 
both of the telluric correction methods and to directly 
compare the results. 

We first compared the standard-star spectrum in each 
IFU with that used by the pipeline routines for the three- 
arm template in each of the spectrographs. Figure 4 
shows the differences between the standard-star spectra 
across the IFUs, where the differences in the Y J-band 
are comparable to those in the i7-band (cf. Fig. 7 from 
Davies, R. et al. 2013). The qualitative agreement be¬ 
tween the IFUs in our region of interest (1.15-1.22/im) 
is generally very good. 

To quantify the difference the two telluric methods 
would make to our analysis, we performed the steps de¬ 
scribed in Section 3.2 for both templates. We then used 
the two sets of reduced science data (reduced with both 
methods of the telluric correction) to compute stellar pa¬ 
rameters for our targets. The results of this comparison 
are detailed in Section 4.1. 

3.2. Telluric Correction Implementation 

To improve the accuracy of the telluric correction, 
for both methods mentioned above, we implemented 
additional recipes beyond those of the KMOS/esorex 
pipeline. These recipes were employed to account for 
two effects which could potentially degrade the quality 
of the telluric correction. The first corrects for any po¬ 
tential shift in wavelength between each science spectrum 
and its associated telluric standard. The most effective 
way to implement this is to cross-correlate each pair of 
science and telluric-standard spectra. Any shift between 
the two is then applied to the telluric standard using a 
cubic-spline interpolation routine. 

The second correction applied is a simple spectral scal¬ 
ing algorithm. This routine corrects for differences in line 
intensity of the most prominent features common to both 
the telluric and the science spectra. To find the optimal 
scaling parameter the following formula is used: 

T 2 = (T 1 +c)/(T 1 -c), (1) 

where T^ is the corrected telluric-standard spectrum, T\ 
is the initial telluric standard spectrum and c is the scal¬ 
ing parameter. 
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Table 1 

Measured velocity resolution and resolving power (. R ) across each detector 


Detector 

IFUs 

Ne A 1.17700 /im 

FWHM (km s -1 ) R 

Ar Al.21430 
FWHM (km s- 1 ) 

/im 

R 

1 

1-8 

88.04 ±2.67 

3408 ± 103 

85.45 ±2.67 

3511 ±110 

2 

9-16 

82.83 ±2.48 

3622±108 

80.30 ±3.05 

3736±142 

3 

17-24 

103.23 ±2.73 

2906 ±77 

101.25 ±2.99 

2963 ±87 


Table 2 

Summary of VLT-KMOS targets in NGC6822 


ID 

S/N 

a (J2000) 

6 (J2000) 

B 

V 

R 

J 

H 

K s 

RV (km s 1 ) 

Notes 

NGC6822-RSG01 

223 

19:44:43.81 

-14:46:10.7 

20.83 

18.59 

17.23 

14.16 

13.37 

13.09 

-69.9 ±3.7 

Sample 

NGC6822-RSG02 

120 

19:44:45.98 

-14:51:02.4 

20.91 

18.96 

17.89 

15.53 

14.72 

14.52 

-66.4 ±6.6 

Sample 

NGC6822-RSG03 

94 

19:44:47.13 

-14:46:27.1 

21.30 

19.41 

18.41 

16.13 

15.35 

15.12 

-43.0 ±6.0 


NGC6822-RSG04 

211 

19:44:47.81 

-14:50:52.5 

20.74 

18.51 

17.22 

14.37 

13.58 

13.30 

-64.5 ±3.1 

LM12 (Ml), Sample 

NGC6822-RSG05 

104 

19:44:50.54 

-14:48:01.6 

20.83 

18.95 

17.97 

15.75 

14.98 

14.79 

—93.4 ± 17.1 


NGC6822-RSG06 

105 

19:44:51.64 

-14:48:58.0 

21.33 

19.45 

18.32 

15.81 

14.95 

14.72 



NGC6822-RSG07 

145 

19:44:53.46 

-14:45:52.6 

20.36 

18.43 

17.38 

15.06 

14.30 

14.08 

-72.6 ±3.5 

LM12 (M4.5), Sample 

NGC6822-RSG08 

103 

19:44:53.46 

-14:45:40.1 

20.88 

19.14 

18.17 

15.95 

15.16 

14.98 

-73.0 ±5.8 

LM12 (K5), Sample 

NGC6822-RSG09 

201 

19:44:54.46 

-14:48:06.2 

20.56 

18.56 

17.35 

14.43 

13.67 

13.34 

-79.9 ±3.7 

LM12 (Ml), Sample 

NGC6822-RSG10 

302 

19:44:54.54 

-14:51:27.1 

19.29 

17.05 

15.86 

13.43 

12.66 

12.42 

-57.2 ±4.7 

LM12 (M0), Sample 

NGC6822-RSG11 

327 

19:44:55.70 

-14:51:55.4 

19.11 

16.91 

15.74 

13.43 

12.70 

12.43 

-67.9 ±3.1 

LM12 (M0), Sample 

NGC6822-RSG12 

100 

19:44:55.93 

-14:47:19.6 

21.43 

19.56 

18.52 

16.14 

15.33 

15.14 

—61.2 ± 6.4 

LM12 (K5) 

NGC6822-RSG13 

106 

19:44:56.86 

-14:48:58.5 

21.05 

19.06 

18.04 

15.81 

15.05 

14.85 

-55.0 ± 10.0 


NGC6822-RSG14 

284 

19:44:57.31 

-14:49:20.2 

19.69 

17.41 

16.20 

13.52 

12.76 

12.52 

— 52.7 ± 4.4 

LM12 (Ml), Sample 

NGC6822-RSG15 

124 

19:44:59.14 

-14:47:23.9 

21.30 

19.17 

18.05 

15.58 

14.74 

14.50 

-74.2 ± 13.4 


NGC6822-RSG16 

107 

19:45:00.24 

-14:47:58.9 

21.27 

19.20 

18.10 

15.60 

14.80 

14.57 

-68.8 ±4.5 


NGC6822-RSG17 

167 

19:45:00.53 

-14:48:26.5 

20.84 

18.75 

17.51 

14.70 

13.86 

13.61 

-62.1 ±4.0 

Sample 

NGC6822-RSG18 

104 

19:45:06.98 

-14:50:31.1 

21.06 

19.12 

18.06 

15.74 

14.94 

14.78 

-87.8 ±9.4 

Sample 


Note. — Optical data from Massey et al. (2007), with typical photometric uncertainty 0.016, 0.006, 0.010 in B,V and R bands, respectively. 
Near-IR data from the UKIRT survey (see Sibbons et al. 2012, for details), with typical errors 0.015, 0.010, 0.012, in J, H and K bands, respectively. 
Targets observed by Levesque & Massey (2012) are indicated by ‘LM12’ in the final column (with their spectral classifications in parentheses). 
Targets used for abundance analysis are indicated by the comment ‘Sample’. 


To determine the required scaling, telluric spectra are 
computed for —0.5 < c < 0.5, in increments of 0.02 
(where a perfect value, i.e. no difference in line strength, 
would be c = 0). Each telluric spectrum is used to correct 
the science data and the standard deviation of the counts 
across the spectral region is computed for each corrected 
spectrum. The minimum value of the standard-deviation 
matrix defines the optimum scaling. For this algorithm, 
only the region of interest for our analysis is considered 
(i.e. 1.15-1.22 /im). 

The final set of telluric-standard spectra from the 
KMOS/esorex reductions were modified using these ad¬ 
ditional routines and were then used to correct the sci¬ 
ence observations for the effects of the Earth’s atmo¬ 
sphere. 

As an alternative to observing telluric standard stars, a 
new telluric correction package, molecfit, allows one to 
calculate a telluric spectrum based on atmospheric mod¬ 
elling. Briefly, the software uses a reference atmospheric 
profile to estimate the true profile for the time and loca¬ 
tion of the science observation. This model is then used 
to create a telluric spectrum which can be used to correct 
the observations. 

This software has been shown to work well on a variety 
of VLT instruments (Smette et al. submitted) and has 
been rigorously tested using X-shooter spectra (Kausch 
et al. submitted). However, the package has yet to be 
tested thoroughly on lower resolution observations such 


as those from KMOS. Our first tests appear encouraging; 
however, pending further characterisation of the KMOS 
data cubes (e.g. small variations in spectral resolving 
power leading to sky residuals, see Section 3.3), we will 
investigate the potential of the molecfit package in fu¬ 
ture papers of this series. 

3.3. Sky Subtraction 

Initial inspection of the extracted stellar spectra re¬ 
vealed minor residuals from the sky subtraction process. 
Reducing these cases with the c sky_tweak’ option within 
the KMOS/esorex reduction pipeline was ineffective to 
improve the subtraction of these features. Any residual 
sky features could potentially influence our results by 
perturbing the continuum placement within the model 
fits, which is an important aspect of the fitting process 
(see Gazak et al. 2014a; Davies et al. 2015, for more dis¬ 
cussion). Thus, pending a more rigorous treatment of the 
data (e.g. to take into account the changing spectral res¬ 
olution across the array), we exclude objects showing sky 
residuals from our analysis. Of the 18 observed targets, 
11 were used to derive stellar parameters (as indicated 
in Table 2). 

3.4. Stellar Radial Velocities 

Radial velocities for each target are listed in Table 2. 
To check the accuracy of the wavelength solution pro¬ 
vided by the data reduction pipeline, this solution is 
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Figure 4. Comparison of J-band spectra of the same standard star in each IFU. The ratio of each spectrum compared with that from 
the IFU used in the three-arm telluric method is shown, with their mean and standard deviation (p, and cr respectively). Red lines indicate 
li = 1.0, a = 0.0 for each ratio. The blue shaded area signifies the region used in our analysis, within which the discrepancies between the 
IFUs are generally small. This is reflected in the standard deviation values when only considering this region. (IFUs 13 and 16 are omitted 
as no data were taken with these IFUs.) 


cross-correlated with the wavelength solution of a spec¬ 
trum of the Earth’s telluric features 1 . In general, this is 
a small correction. The science spectra are then telluric 
corrected in the manner described above. 

Radial velocities are estimated by measuring the peak 
of the cross-correlation function between the telluric- 
corrected science spectra (taking the spectrum from the 
brightest spatial pixel) and an appropriate synthetic 
RSG spectrum. The errors are calculated by taking the 
dispersion of each line used to derive the radial velocities 
scaled by the number of lines used, following Lapenna et 
al. (2015). A robust radial velocity could not be derived 
for one of our candidates (NGC6822-RSG06) owing to 
strong residual sky features; therefore, it is not consid¬ 
ered further. 

Radial velocity estimates are shown as a function of 
distance to the centre of NGC6822 in Figure 5. The 
average radial velocity for our targets is —68 =b 12 km s _1 , 
in good agreement with the systemic radial velocity of the 
HI disk (—57±2kms _1 ; Koribalski et al. 2004). Our 

1 Retrieved from http://www.eso.org/sci/facilities/paranal/ de¬ 
commissioned / isaac / tools / spectroscopic_standards. html 


radial velocities also agree with estimates for the two A- 
type supergiants from Venn et al. (2001). This result 
confirms that our candidates are NGC6822 members. 

4. RESULTS 

Stellar parameters (metallicity, effective temperature, 
surface gravity and microturbulence) have been derived 
using the J-band analysis technique described by Davies 
et al. (2010) and demonstrated by Gazak et al. (2014a) 
and Davies et al. (2015). To estimate physical param¬ 
eters this technique uses a grid of synthetic spectra 
to fit observational data, in which the models are de¬ 
graded to the resolution of the observed spectra (Ta¬ 
ble 1). Model atmospheres were generated using the 
MARCS code (Gustafsson et al. 2008) where the range 
of parameters is defined in Table 3. The precision of the 
models is increased by including departures from LTE in 
some of the strongest Fe, Ti and Si atomic lines (Berge- 
mann et al. 2012, 2013). The two strong magnesium lines 
in our diagnostic spectral region are excluded from the 
analysis at present as these lines are known to be affected 
strongly by non-LTE effects (see Figure 6, where the two 
Mg I lines are systematically under- and over-estimated). 
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R (arcseconds) 



Figure 5. Radial velocities of targets shown against their dis¬ 
tance from the galaxy centre. The average radial velocity for the 
sample is —68 d= 12 km s —1 . The green dashed line indicates the 
HI systemic velocity (—57±2km s -1 ; Koribalski et al. 2004). The 
radial velocities of two A-type supergiants from Venn et al. (2001) 
are shown in blue. 


Table 3 

Model grid used for analysis 


Model Parameter 

Min. 

Max. 

Step size 

T ef f (K) 

3400 

4000 

100 


4000 

4400 

200 

[Z] (dex) 

-1.50 

1.00 

0.25 

log 9 (cgs) 

-1.0 

1.0 

0.5 

£ (km s -1 ) 

1.0 

6.0 

1.0 


The non-LTE effects on the formation of the Mg I lines 
will be explored by Bergemann et al. (submitted). 

4.1. Telluric Comparison 

We used these Science Verification data to determine 
which of the two telluric standard methods is most ap¬ 
propriate for our analysis. Table 4 details the stellar 
parameters derived for each target using both telluric 
methods, and these parameters are compared in Figure 7. 
The mean difference in metallicity from the two methods 
is A[Z] = 0.04 ±0.07. Therefore, for our analysis, there 
is no significant difference between the two telluric ap¬ 
proaches. 

4.2. Stellar Parameters and Metallicity 

Table 4 summarises the derived stellar parameters. For 
the remainder of this paper, when discussing stellar pa¬ 
rameters, we adopt those derived using the 24-arm tel¬ 
luric method (i.e. the results in the left-hand part of 
Table 4). The average metallicity for our sample of 11 
RSGs in NGC6822 is [Z] = -0.52 ±0.21. This result is 
in good agreement with the average metallicity derived in 
NGC 6822 from blue supergiant stars (BSGs; Muschielok 
et al. 1999; Venn et al. 2001). 

A direct comparison with metallicities from BSGs is le¬ 
gitimate as the results derived here yield a global metal¬ 
licity ([Z]) which closely resembles the Fe/H ratio derived 
by Venn et al. (2001). While our [Z] measurements are 


also affected by Si and Ti, we assume [Z] = [Fe/H] for 
the purposes of our discussion. Likewise, we can compare 
oxygen abundances (relative to solar) obtained from HII 
regions as a proxy for [Z] by introducing the solar oxygen 
abundance 12 ±log (O/H)© = 8.69 (Asplund et al. 2009) 
through the relation [Z] = 12 ± log (O/H) — 8.69. 

The RSG and BSG stages are different evolutionary 
phases within the life cycle of a massive star, while HII re¬ 
gions are the birth clouds which give rise to the youngest 
stellar population. As the lifetimes of RSGs and BSGs 
are < 50Myr, their metallicity estimates are also ex¬ 
pected to be representative of their birth clouds. 

To investigate the spatial distribution of chemical 
abundances in NGC 6822, in Figure 8 we show the metal¬ 
licities of our RSGs as a function of radial distance from 
the centre of the galaxy, as well as the results from 
Venn et al. (2001) and the indicative estimates from 
Muschielok et al. (1999). 

A least-squares fit to the KMOS results reveals a 
low-significance abundance gradient within the central 
1 kpc of NGC 6822 of —0.5 ± 0.4 dex kpc -1 . The extrap¬ 
olated central metallicity from the fit (i.e. at R = 0) of 
[Z] = —0.30 ±0.15 derived remains consistent with the 
average metallicity assuming no gradient. 

Figure 9 shows the location of our sample in the 
Hertzsprung-Russell (H-R) diagram. Bolometric correc¬ 
tions were computed using the calibration from Davies, 
B. et al. (2013) to calculate luminosities. This figure 
shows that the temperatures derived using the J-band 
method are systematically cooler than the end of the 
evolutionary models (which terminate at the end of the 
carbon-burning phase for massive stars) for Z = 0.002 
from Georgy et al. (2013). This is discussed in Sec¬ 
tion 5.2. 

5. DISCUSSION 

5.1. Metallicity Measurements 

_We find an average metallicity for our sample of 
[Z] = —0.52 ±0.21 which agrees well with the results de¬ 
rived from BSGs (Muschielok et al. 1999; Venn et al. 
2001; Przybilla 2002) and HII regions (Lee et al. 2006). 

We also find evidence for a low-significance metal¬ 
licity gradient within the central lkpc of NGC 6822 
(—0.5 ± 0.4 dex kpc -1 with a x^ed = 1*16; see Figure 8). 
The gradient derived is consistent with the trend re¬ 
ported in Venn et al. (2001) from their results for the two 
BSGs compared with HII regions from Pagel et al. (1980) 
and two planetary nebulae from Richer & McCall (1995) 
at larger galactocentric distances. Our result is also con¬ 
sistent with the gradient derived from a sample of 49 local 
star-forming galaxies (Ho et al. submitted). Including 
the results for BSGs from Venn et al. (2001) in our anal¬ 
ysis gives a consistent gradient (—0.48 ± 0.33 dex kpc -1 ) 
with a Xred = 1-06. Results from Muschielok et al. (1999) 
are not included in the fit as these measurements were 
qualitative estimates of metallicity. 

In contrast, Lee et al. (2006) used the oxygen abun¬ 
dances from 19 HII regions and found no clear evidence 
for a metallicity gradient. Using a subset of the highest 
quality HII region data available, these authors found a 
gradient of —0.16 ± 0.05 dexkpc -1 . Including these re¬ 
sults into our analysis degrades the fit and changes the 
derived gradient significantly (—0.18±0.05dexkpc -1 ; 
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Figure 6. KMOS spectra of the NGC6822 RSGs and their associated best-fit model spectra (black and red lines, respectively). The 
lines used for the analysis from left to right by species are: Fe IAA1.188285, 1.197305, Si IAA1.198419, 1.199157, 1.203151, 1.210353, 
TiIAA1.189289, 1.194954. The two strong Mgl lines are also labelled, but are not used in the fits (see Section 4). 
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Figure 7. Comparison of the model parameters using the two different telluric methods. In each panel, the x-axis represents stellar 
parameters derived using the three-arm telluric method and the y-axis represents those derived using the 24-arm telluric method. Top left: 
metallicity ([Z]), mean difference < A[Z] > = 0.04 ±0.07. Top right: effective temperature (T e ff), mean difference < AT e ff > = —14 ±42. 
Bottom left: surface gravity (log ( 7 ), mean difference < A log (7 > = —0.06 ±0.12. Bottom right: microturbulence (£), mean difference 
< Ad; > = —0.1 ±0.1. In all cases, the distributions are statistically consistent with a one-to-one ratio (black lines). 


Table 4 

Fit parameters for reductions using the two different telluric methods. 


Target 

IFU 


24-arm 

Telluric 


three- 

-arm Telluric 


Teff (K) 

log 3 

£ (km s _i ) 

[z] 

Teff (K) 

log 3 

£ (km s i ) 

[Z] 

NGC6822-RSG01 

6 

3790 ± 80 

-0.0 ±0.3 

3.5 ±0.4 

-0.55 ±0.18 

3860 ± 90 

—0.1 ± 0.5 

3.5 ±0.4 

-0.61 ±0.21 

NGC6822-RSG02 

11 

3850 ±100 

0.4 ± 0.5 

3.5 ±0.4 

-0.78 ±0.22 

3810 ±110 

0.4 ± 0.5 

3.3 ± 0.5 

-0.65 ±0.24 

NGC6822-RSG04 

12 

3880 ± 70 

0.0 ±0.3 

4.0 ±0.4 

-0.32 ±0.16 

3880 ± 70 

0.0 ±0.3 

4.0 ±0.4 

-0.32 ±0.16 

NGC6822-RSG07 

2 

3970 ± 60 

0.4 ± 0.5 

3.9 ±0.4 

-0.58 ±0.19 

3990 ± 80 

0.1 ± 0.5 

3.8 ± 0.5 

-0.56 ±0.14 

NGC6822-RSG08 

3 

3910 ±100 

0.6 ± 0.5 

3.0 ±0.4 

-0.58 ±0.24 

3910 ±100 

0.6 ± 0.5 

3.0 ±0.4 

-0.58 ±0.24 

NGC6822-RSG09 

4 

3980 ±60 

0.1 ± 0.4 

3.7 ± 0.4 

-0.38 ±0.16 

3990 ± 80 

—0.1 ± 0.5 

3.6 ±0.4 

—0.44 ± 0.17 

NGC6822-RSG10 

14 

3900 ±80 

—0.3 ± 0.5 

3.7 ± 0.4 

—0.67 ± 0.16 

3850 ± 80 

—0.3 ± 0.5 

3.5 ±0.4 

-0.59 ±0.19 

NGC6822-RSG11 

15 

3870 ± 80 

—0.4 ± 0.5 

4.2 ±0.5 

-0.53 ±0.19 

3850 ± 60 

-0.3 ±0.4 

4.3 ±0.5 

-0.49 ±0.17 

NGC6822-RSG14 

17 

3910±110 

-0.5 ±0.5 

3.6 ±0.5 

—0.20 ± 0.21 

3880±110 

-0.5 ±0.5 

3.6 ±0.5 

-0.15 ±0.24 

NGC6822-RSG17 

21 

3890±120 

0.1 ± 0.5 

3.0 ±0.4 

-0.43 ±0.28 

3890±120 

0.1 ± 0.5 

3.0 ±0.4 

-0.43 ±0.28 

NGC6822-RSG18 

18 

3810±130 

0.4 ± 0.5 

2.2 ± 0.4 

-0.68 ±0.31 

3740 ±130 

0.4 ± 0.5 

2.1 ± 0.5 

-0.54 ±0.41 
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Figure 8. Derived metallicities for 11 RSGs in NGC6822 shown 
against their distance from the galaxy centre; the average metal- 
licity is [Z] = —0.52 ± 0 . 21 . Blue and red points show BSG re¬ 
sults from Venn et al. (2001) and Muschielok et al. (1999), re¬ 
spectively. A least-squares fit to the KMOS results reveals a low- 
significance abundance gradient (see text for details). For compar¬ 
ison, R 25 = 460" (= 1.03kpc; McConnachie 2012 ). 
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Figure 9. H-R diagram for the 11 RSGs in NGC6822. Evolu¬ 
tionary tracks including rotation (v/v c = 0.4) for SMC-like metal- 
licity (Z = 0 . 002 ) are shown in grey, along with their zero-age mass 
(Georgy et al. 2013). Bolometric corrections are computed using 
the calibration from Davies, B. et al. (2013). We note that, com¬ 
pared with the evolutionary tracks, the observed temperatures of 
NGC6822 RSGs are systematically cooler. This is discussed in 
Section 5.2. 


Xred = 1-78). At this point it is not clear whether the in¬ 
dication of a gradient obtained from the RSGs and BSGs 
is an artefact of the small sample size, or indicates a dif¬ 
ference with respect to the HII region study. 

There have been a number of studies of the metallic- 
ity of the older stellar population in NGC6822. From 
spectroscopy of red giant branch (RGB) stars, Tolstoy et 
al. (2001) found a mean metallicity of [Fe/H] = —0.9 
with a reasonably large spread (see their Fig. 19). 
More recently, Sibbons et al. (2012) derived metallic- 
ities using a population of AGB stars within the cen¬ 
tral 4 kpc of NGC 6822. They found an average metal¬ 
licity of [Fe/H] = —1.29 ± 0.07dex. Likewise, Kirby 


et al. (2013) used spectra of red giant stars within 
the central 2 kpc and found an average metallicity of 
[Fe/H] = —1.05 ±0.49 dex. We note that none of these 
studies found any compelling evidence for spatial varia¬ 
tions in the stellar metallicities, which is not surprising 
given that, in disc galaxies, radial migration is thought 
to smooth out any abundance gradients over time. The 
stellar populations used for these studies are known to be 
substantially older than our sample; therefore, owing to 
the chemical evolution in the time since the birth of the 
older populations, we expect the measured metallicities 
to be significantly lower. 

The low metallicity of the young stellar population and 
the interstellar medium (ISM) in NGC 6822 can be un¬ 
derstood as a consequence of the fact that it is a relatively 
gas-rich galaxy with a mass of M hi = 1.45 x 10 8 M© 
(Koribalski et al. 2004) and a total stellar mass of 
M* = 0.83 - 1.70 x 10 8 M© (Woo et al. 2008; Kirby et 
al. 2013; Weisz et al. 2014). 

The simple closed-box chemical-evolution model re¬ 
lates the metallicity mass fraction Z(t ) at any time to 
the ratio of stellar to gas mass through: 


Z(t) 


y 

1 -R 


In 


M* (t )" 

M g (t)\ ’ 


( 2 ) 


where y is the fraction of metals per stellar mass 
produced through stellar nucleosynthesis (the so-called 
yield) and R is the fraction of stellar mass returned to 
the ISM through stellar mass-loss. 

According to Kudritzki et al. (in prep.), the ratio 
y/( 1 — R ) can be empirically determined from the fact 
that the metallicity of the young stellar population in 
the solar neighbourhood is solar, with a mass fraction of 
Z© = 0.014 (Nieva & Przybilla 2012). With a stellar-to¬ 
gas mass column density of 4.48 in the solar neighbour¬ 
hood (Wolfire et al. 2003; Bovy & Rix 2013) one then 
obtains y/( 1 — R) = 0.0082 = 0.59Z© with an uncer¬ 
tainty of 15% dominated by the 0.05 dex uncertainty of 
the metallicity determination of the young population. 

Accounting for the presence of helium and metals in 
the neutral interstellar gas we can turn the observed HI 
mass in NGC 6822 into a gas mass via M^ = 1.36 M hi 
and use the simple closed-box model to predict a metal¬ 
licity of [Z] = —0.44 to —0.69, in good agreement with 
our value obtained from RSG spectroscopy. 

As discussed above, the older stellar population of 
AGB stars has a metallicity roughly 0.8 dex lower than 
the RSGs. In the framework of the simple closed-box 
model this would correspond to a period in time where 
the ratio of stellar to gas mass was ~0.1 (much lower 
than the present value of 0.42 to 0.86). The present star- 
formation rate of NGC 6822 is ^0.02 M©yr -1 (Gratier 
et al. 2010; Efremova et al. 2011). At such a level of 
star formation it would have taken 5 Gyr to produce the 
presently observed stellar mass and to arrive at the av¬ 
erage metallicity of the young stellar population from 
that of the AGB stars (of course, again relying on the 
simple closed-box model). Evidence suggests that the 
star-formation rate was substantially lower in the past 
(Efremova et al. 2011; Weisz et al. 2014), therefore, the 
build up of the observed stellar mass would have taken 
correspondingly longer. 
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Figure 10. Effective temperatures as a function of metallicity 
for four different data sets using the J-band analysis technique. 
There appears to be no significant variation in the temperatures 
of RSGs over a range of 0.55 dex. These results are compiled from 
the LMC, SMC (blue and red points, respectively; Davies et al. 
2015), Perseus OB-1 (a Galactic RSG cluster; cyan; Gazak et al. 
2014a) and NGC6822 (green). Mean values for each data set are 
shown as enlarged points in the same style and colour. The x-axis 
is reversed for comparison with Figure 11. 

Given the irregularities present in the stellar and 
gaseous morphology of NGC6822, this galaxy may not 
be a good example of a closed-box system; however, it 
is remarkable that the closed-box model reproduces the 
observed metallicity so closely. 

5.2. Temperatures of RSGs 

Effective temperatures have been derived for 11 RSGs 
from our observed sample in NGC6822. To date, this 
represents the fourth data set used to derive stellar pa¬ 
rameters in this way and the first with KMOS. The previ¬ 
ous three data sets which have been analysed are those of 
11 RSGs in Perseus OB-1, a Galactic star cluster (Gazak 
et al. 2014a), nine RSGs in the LMC and 10 RSGs in the 
SMC (both from Davies et al. 2015). These results range 
from Z = Z 0 in Perseus OB-1 to Z = 0.3 Z 0 in the SMC, 
around 0.5 dex in metallicity. 

We compare the effective temperatures derived in this 
study with those of the previous results in different en¬ 
vironments. Their distribution is shown as a function of 
metallicity in Figure 10. Additionally, Figure 11 shows 
the H-R diagram for the four sets of results. Bolometric 
corrections to calculate the luminosities for each sample 
were computed using the calibration from Davies, B. et 
al. (2013). 

From these figures, we see no evidence for significant 
variations in the average temperatures of RSGs with re¬ 
spect to metallicity. This is in contrast with current evo¬ 
lutionary models which display a change of ^450 K, for 
a M = 15M 0 model, over a range of solar to SMC-like 
metallicities (Ekstrom et al. 2012; Georgy et al. 2013). 

For solar metallicity, observations in Perseus OB-1 are 
in good agreement with the models (see Fig. 9 in Gazak 
et al. 2014a). However, at SMC-like metallicity, the 
end-points of the models are systematically warmer than 
the observations. The temperature of the end-points of 
the evolutionary models of massive stars could depend 



Figure 11. H-R diagram for RSGs in Perseus OB-1 (cyan), LMC 
(blue), SMC (red) and NGC 6822 (green) which have stellar param¬ 
eters obtained using the J-band method. This figure shows that 
there appears to be no significant temperature differences between 
the four samples. Solid grey lines show SMC-like metallicity evo¬ 
lutionary models including rotation (Georgy et al. 2013). Dashed 
grey lines show solar metallicity evolutionary models including ro¬ 
tation (Ekstrom et al. 2012). 

on the choice of the convective mixing-length parameter 
(Schaller et al. 1992). The fact that the models produce 
a higher temperature than observed could imply that the 
choice of a solar-like mixing-length parameter does not 
hold for higher mass stars at lower metallicity. 

Lastly, we note that the average spectral type of RSGs 
tends towards an earlier spectral type with decreasing 
metallicity over this range (Humphreys 1979; Levesque 
& Massey 2012). We argue that this is not in contradic¬ 
tion to the above results. Spectral types are derived for 
RSGs using the optical TiO band-heads at ^0.65 /im, 
whereas in this study temperatures are derived using 
near-IR atomic features (as well as the line-free pseudo¬ 
continuum). The strengths of TiO bands are dependent 
upon metallicity which means that the spectral classi¬ 
fication for RSGs at a constant temperature will differ 
(Davies, B. et al. 2013). Therefore, although historically 
spectral type has been used as a proxy for temperature, 
this assumption does not provide an accurate picture for 
RSGs in differing environments. 

5.3. AGB Contamination 

As mentioned in Section 2.1, massive AGB stars are 
potential contaminants to our sample. These stars have 
similar properties to RSGs and can occupy similar mass 
ranges as lower-mass RSGs (Herwig 2005); however, their 
lifetimes are around ^250 Myr (Doherty et al. 2010). 
Wood et al. (1983) argued for an AGB luminosity limit 
(owing to the limit on the mass of the degenerate core) of 
Mboi ~ —7.1. Using this maximum luminosity, corrected 
for the distance to NGC 6822, yields K = 14.0. Four of 
our analysed stars have iGband magnitudes fainter than 
this limit, but excluding the results for these does not 
significantly alter any of our results. However, given the 
near-IR luminosity and colours of these stars, if indeed 
they are AGB stars, they are probably younger massive 
AGB stars (e.g. Siess 2008) and therefore would arguably 
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still trace the relatively young stellar population. 

6 . CONCLUSIONS 

KMOS spectroscopy of RSG stars in NGC 6822 is pre¬ 
sented. The data were telluric corrected in two differ¬ 
ent ways and the standard three-arm telluric method is 
shown to work as effectively (in most cases) as the more 
time expensive 24-arm telluric method. Radial velocities 
of the targets are derived and are shown to be consis¬ 
tent with previous results in NGC 6822, confirming their 
extragalactic nature. 

Stellar parameters are calculated for 11 RSGs using 
the J-band analysis method outlined by Davies et al. 
(2010). The average metallicity within NGC6822 is 
[Z] = —0.52 ±0.21, consistent with previous abundance 
studies of young stars. We find an indication of a metal¬ 
licity gradient within the central 1 kpc, but with a low 
significance caused by the small size and limited spatial 
extent of our RSG sample. To conclusively assess the 
presence of a metallicity gradient among the young pop¬ 
ulation within NGC 6822 a larger systematic study of 
RSGs is needed. 

The chemical abundances of the young and old stellar 
populations of NGC 6822 are well explained by a simple 
closed-box chemical evolution model. However, while an 
interesting result, we note that the closed-box model is 
unlikely to be a good assumption for this galaxy given 
its morphology. 

The effective temperatures of those in this study are 
compared with RSGs analysed using the same method 
in previous studies. Using results which span 0.55 dex 
in metallicity (solar to SMC) within four galaxies, we 
find no evidence for a systematic variation in average 
effective temperature with respect to metallicity. This is 
in contrast with evolutionary models which, for a similar 
change in metallicity, predict a shift in the temperature 
of RSGs of up to 450 K. We argue that an observed shift 
in average spectral type of RSGs over this metallicity 
range does not imply a shift in average temperature. 

These observations were taken as part of the KMOS 
Science Verification programme. With guaranteed time 
observations we have obtained data for RSGs in NGC 300 
and NGC 55 at distances of ^1.9 Mpc, as well as obser¬ 
vations of super-star clusters in M 83 and the Antennae 
galaxy at 4.5 and 20Mpc, respectively. Owing to the 
fact that RSGs dominate the light output from super- 
star clusters (Gazak et al. 2013) these clusters can be 
analysed in a similar manner (Gazak et al. 2014b), which 
will provide metallicity measurements at distances a fac¬ 
tor of 10 larger than using individual RSGs! This work 
is the first step towards an ambitious proposal to survey 
a large number of galaxies in the Local Volume, moti¬ 
vated by the twin goals of investigating their abundance 
patterns, while also calibrating the relationship between 
galaxy mass and metallicity in the Local Group. 

We thank the referee for a thorough review and con¬ 
structive comments. We thank Mike Irwin for provid¬ 
ing the photometric catalogue from the WFCAM ob¬ 
servations. RPK and JZG acknowledge support by the 
National Science Foundation under grant AST-1108906. 
The NSO/Kitt Peak FTS high-resolution telluric spec¬ 
tra used here to investigate the wavelength solution were 


produced by NSF/NOAO. 
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